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An analytical method is presented to solve generalized QCD evolution equa- 
tions for the time development of parton cascades in a nuclear environment. 



In addition to the usual parton branching processes in vacuum, these evolu- 
tion equations provide a consistent description of interactions with the nuclear 



Qh. medium by accounting for stimulated branching processes, fusion and scat- 



tering processes that are specific to QCD in medium. Closed solutions for the 
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momentum and virtuality are obtained under some idealizing assumptions 
about the composition of the nuclear medium. Several characteristic features 
of the resulting parton distributions are discussed. One of the main conclu- 
sions is that the evolution of a parton shower in medium is dilated as compared 
to free space and is accompanied by an enhancement of particle production. 
These effects become stronger with increasing nuclear density. 

PACS Indices: 25.75.+r, 12.38.Bx, 12.38.Mh, 24.85. +p 

e-mail: klaus@suryall.cern.ch 

CERN-TH. 7206/94, March 1994 





1. INTRODUCTION 

In a preceding paper |T[, B. Miiller and myself addressed the question, how is "QCD in 
medium" modified as compared to "QCD in vacuum" when one attempts to extend well 
developed perturbative QCD techniques from high energy hadron-hadron collisions to 
ultra-relativistic heavy ion collisions. The investigation of [Q] condensed in the derivation of 
a generalized form of QCD evolution equations that describe the time evolution of parton 
distributions in a nuclear medium, solutions of which I will present in this paper. This 
aspect is of great importance for the future ultra-relativistic heavy ion collider experiments 
at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider 
(LHC), where new phenomena associated with "QCD in medium" are expected to clearly 
modify naive extrapolations from hadron-hadron collisions. Although in the last couple 
of years, numerical simulations with QCD based Monte Carlo models have provided 
considerable new insight in the microscopic parton dynamics of ultra-relativistic heavy ion 
collisions, one is still far from a complete picture and a truly quantitative description. Aside 
from the rather intransparent complexity of these computer simulations, one still relies on 
a substantial amount of phenomenological modeling of nuclear and dense medium effects 
due the current lack of better knowledge about the details of such mechanisms. This is 
reflected by a very large uncertainty in quantitative predictions: although the different 
model calculations generally agree very well for pp collisions at collider energies, they differ 
in their predictions, for e.g. charged particle multiplicities, in heavy ion AA collisions by a 
factor of 2 or more. 

These issues lead to the main motivation for this work, namely the necessity to go back 
to a simpler, but more fundamental and transparent level, and study within the firmly es- 
tablished framework of perturbative QCD (rather than extending phenomenological model 
input) the microscopic dynamics of quarks and gluons in the hot, ultra-dense environment 
that may be created in heavy ion collisions |7|-|J]. In order to do so, it is inevitable to dis- 
entangle the various nuclear and medium effects from each other by investigating physical 



1 



situations where such medium phenomena may be singled out, both in theoretical and exper- 
imental research. For instance, deep inelastic scattering on heavy nuclei ||10||, charmonium 



production by Drell-Yan processes in nuclei [11|, energy loss of partons, jet quenching and 



multiple scattering phenomena in hadron- nucleus collisions ||12|| , etc., provide opportunities 
to explore nuclear modifications to the QCD parton picture that has been so successful in 
high energy particle physics. 

As a first step in this direction, in Ref. a generalized form of the well known Dokshitzer- 



Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations [|13| |15|1 was derived. In the 
framework of the leading logarithmic approximation (LLA) |T6|-[19|, we obtained a set of 
coupled integro-differential equations for the time evolution (rather than the Q 2 -evolution) 
of the off-shell parton distribution functions in nuclear matter. In addition to the usual 
parton branching processes in vacuum, these evolution equations provide a self-contained 
description of interactions with the nuclear medium by accounting for stimulated branching 
processes, fusion and scattering processes, which are specific to QCD in medium. This ap- 
proach resulted in a probabilistic description within perturbative QCD of the time variation 
of parton distributions as in non-equilibrium kinetic theory, however now including consis- 
tent treatment of off-shell propagation of partons, a concept that is absent in semiclassical 
transport theory [^0|. The essential condition for the applicability of perturbative techniques 



is that the QCD factorization theorem J2]J holds also in the presence of dense nuclear mat- 
ter, in which multiple interactions of the evolving partons with this nuclear medium may 
be separated in space-time by only very short distances. Recent investigations by McLerran 



and Venugopalan |22| are encouraging in this direction, since they showed that a consistent 
perturbative calculation of parton structure functions at small values of the Bjorken variable 
x becomes possible when one considers the limit of a very thick nuclear target. The con- 
dition for the applicability of perturbative QCD is then that medium induced effects, such 
as color screening and rescattering, provide dynamical cut-offs on a scale short compared to 
the QCD renormalization scale A . 

Exploring this insight, the present paper as a sequel to the work of Ref. [Q employs 



the probabilistic parton cascade picture [|15"1 , [2"3| to describe the evolution of a parton shower 
inside nuclear matter, triggered by a highly excited initial quark or gluon. Here I will start 
from the generalized evolution equations and solve for the parton spectra with respect to the 
variables time, longitudinal momentum and virtuality, under some idealizing assumptions 
about the composition of the nuclear medium. Before going to the heart of the matter, 
let me remark the following: The approach is at this point in a yet idealized stage, since 
various unsettled problems are left out, and a number of simplifying approximations are 
made. Still, it is a starting point which can be step by step improved, similarly maybe as 
the development of perturbative QCD in hard processes that started some 15 years ago. The 
main assumptions and approximations in Ref. |]] and in the present paper can be stated as 
follows: 

(i) Factorization of short distance interactions of shower partons with the nuclear medium 

from long range non-perturbative forces is assumed. This means, even in the presence 
of dense medium where a parton can encounter multiple successive interactions, a 
probabilistic description of local, non-interfering interactions applies at sufficiently 
high energies. 

(ii) Neglect of interference effects that affect particularly the softest partons: destructive 

interference of interaction amplitudes in both, coherent succesive small angle emis- 
sions ('angular ordering'), and multiple sequential scatterings ('Landau- Pomerantchuk- 
Migdal effect'). 

(iii) Nuclear effects due to long range correlations, e.g. nuclear shadowing (antishadowing) 
not associated with truly perturbative parton interactions are ignored. 

The remainder of the paper is organized as follows. Sec. 2 is devoted to recall the intuitive 
picture developed in Ref. |T[j of the probabilistic parton evolution in medium. The theoretical 
framework is reviewed and the master equations describing the time development of the 
distribution of off-shell partons in a nuclear environment are summarized. Sec. 3 then deals 
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with linearizing and solving these equations analytically by making a number of simplifying 
assumptions regarding the nuclear medium. Explicit closed solutions for parton showers 
initiated by an energetic time-like quark or gluon are obtained, which in the special case of 
absence of nuclear matter reduces to the well known LLA solutions for parton distributions, 
however now time being the evolution variable rather than virtuality. Some general features 
concerning the x-dependence as well as the influence of the nuclear density are discussed. A 
summary and outlook is given in Sec. 4. 

2. THE COUPLED EVOLUTION EQUATIONS FOR QUARKS, ANTI- 
QUARKS AND GLUONS IN MEDIUM 

Let me start by recalling the essential considerations of Ref. [|J that lead to the derivation 
of a coupled set of integro-differential equations for a system of off-shell quarks (antiquarks) 
Qi (ft) with flavors i = l,...,nf, and gluons g. As will become clear, these Boltzmann 
type equations describe the time evolution of parton cascades in a nuclear medium with 
account of off-shell propagation of time-like virtual partons. Although as mentioned before, 
the evolution equations may be applied to a variety of physical situations, let me specify 
here for the purpose of lucidity a definite, though somewhat idealized scenario: Consider a 
fast parton injected into infinitely extended nuclear matter by some highly localized process 
of space- time extent, {Ql)~ 112 < A -1 . This primary parton then propagates through the 
nucleus and initiates a cascade or shower of secondary partons. The partons in the cascade 
can either radiate bremsstrahlung gluons or produce gg-pairs, collide with partons of the 
nucleus (the medium), or absorb nuclear partons. I will call such an event a parton shower 
or parton cascade. I emphasize that I will explicitly distinguish between the shower partons 
on the one hand, and the nuclear partons on the other hand which initially are coherently 
bound in the wave function of the nucleus. Also, for the time being, I consider here the 
evolution of a single cascade and therefore only account for interactions of the cascade with 
the medium, but neglect possible interactions between simultaneously evolving cascades. 

The prototype reaction of such a scenario is a proton-nucleus (pA) collision at very 
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high beam energy, where one may trigger on a high energy time-like quark or gluon that is 
produced at some point of time to inside the heavy nucleus by a hard scattering, and single 
out the development of the resulting parton shower inside the nucleus with a sufficiently 
homogenous spatial density. For the description of such a parton cascade it is convenient 
to choose the nucleon-nucleon center- of -mass frame (CMjvat) [23] in which each nucleon has 
the same value of longitudinal momentum P, 

P^ = +P , Pj A ) = -AP , p[ p) = p[ A) = , (1) 

with 



= ^AAP 2 + M 2 N (1 + A 2 ) , = 2P , (2) 

where A is the nuclear mass number, the nucleon mass, and P/M^ ^> 1 is assumed. 
For example, at the RHIC (LHC) the maximum available beam energy implies P/Mn > 
100 (3000) even for p + Au collisions. I will describe the longitudinal evolution of a parton 
shower along the shower axis (z-axis), which I define parallel to the direction of momentum 
of the initiating primary parton. It is convenient to parametrize the four-momenta k = k^ = 
(E, k z , k±) of the shower partons such that for the primary parton, 

k = (*o p +l^p-> x p , o) , (3) 

whereas for the j th secondary parton, 

/ O 2 + k 2 \ 

kj = \x j P+^-^-; zjP, k. j , (4) 

where Xj = p Z j/P can be either positive or negative depending on the partons direction 
of propagation along the z-axis. Furthermore, k\- < Q 2 P 2 is assumed, and all rest 
masses are neglected. It is important to realize that energy and momentum are independent 
variables, since one is dealing with off-shell particles of virtuality Q 2 > with a continuous 
invariant mass distribution. The evolution of the parton system is then described by the 
change of parton number densities, which are defined as 



a(x,k 2 ± ;Q 2 ,t) = jf dQ' 2 d ^J^Q l2 > ( a = ?*» # 5 * = 1, . . . , n f ) , (5) 

or, when integrated over transverse momentum, 

a(x,Q 2 ,i) = J dk 2 x a(x,k 2 ± ;Q 2 ,t) . (6) 

The relation between the parton number densities fl5|) to the single-particle phase-space 
distributions is expressed by the number of partons of type a present at time t 

N a {t) = Jd 3 rJd 4 kF a (E,k;t,r) = J d A k f a (E,k;t) , (7) 

with F a denoting the corresponding phase-space density of off-shell partons with virtuality 
Q 2 = E 2 — k 2 in the phase-space volume dEd 3 kd 3 r around k^ and r at time t, and f a in the 
second line representing the spatially integrated energy-momentum distribution. 

In order to describe the time evolution of the parton densities ([5]) or (Q) one has to relate 
the change of the dynamical variables x and Q 2 with the laboratory time t which plays the 
role of an external 'parameter' rather than being an intrinsic kinematical quantity. In Ref. 
U it was shown, using time-dependent perturbation theory, that for the case of a parton 
shower evolving in vacuum (i.e. by successive branchings only), one finds that the time scale 
for a branching chain, starting at to = with initial values Ql and xo, is in the average given 
by 

with (zj = Xj/xj-i) 



x = x n = zxz-i . . . z n x , Q 2 = Q 2 n < Q 2 n _ x < . . . < Ql , (9) 

in accord with the uncertainty principle. Thus, in the absence of interactions with surround- 
ing matter, the time variable t is related to x and and the degree of off-shellness Q 2 by the 
typical life-time 1/Q in the parton's restframe, boosted by the Lorentz factor 7 = \x\P/Q. 
However, when considering a parton shower inside nuclear matter, scattering and fusion 
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processes with nuclear partons will compete with spontaneous branchings and consequently 
disrupt a steady decrease of average virtuality with time. Both scatterings and fusions can 
increase the virtuality of a parton by energy-momentum transfer. This means that the re- 
lation (|D between x, Q 2 and t cannot be deterministic anymore, because interactions with 
the nuclear medium occur stochastically according to the density of nuclear partons. Each 
scattering or fusion may "rejuvenate" a shower parton by providing it with a virtuality 
Q' 2 > Q 2 so that the parton after the interaction appears to be younger (t < t'). To keep 
track of these repeated rejuvenations, one has to reset the clock for the particular shower 
parton after each such interaction with the medium. For this purpose a new independent 
dynamic variable r, called the age of a parton, was introduced in Ref. QXJ] . Instead of @, 
one finds 

r = r(x,t) = ^jpj , (10) 

which determines the typical age of a parton depending on x, Q 2 and implicitely also on 
t, and which introduces an additional time scale that reflects the external influence of the 
medium on the time evolution of the parton cascade. Correspondingly, the parton number 
densities (|5]) and (^) must be generalized as 

.(*,r, t )=jT ir'jikl^,, (11) 

which now has the meaning of the probability density for finding a parton of type a at time 
t with age r = \x\P/Q 2 and momentum fraction x. The age variable measures here the 
influence of the nuclear medium on the development of the parton cascade. Note that in 
absence of a background medium the age variable looses its independent character, i.e. it 
then evolves parallel to laboratory time and r oc t. I will return to that point in Sec. 3. 

2.1 The evolution equations 

For the remainder of this paper I will be concerned only with transverse momentum inte- 
grated quantities, so that it suffices to describe the evolution in terms of the parton number 
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densities g«, (<&) and g, denned by eq. ([Ll]). In addition, it is convenient to introduce the 
parton momentum densities Qi (Qi) and G 



Qi(x,r,t) = xqi(x,r,t) 



G(x,r,t) = xg(x,r,t) 



(12) 



i.e. the parton number densities weighted with the longitudinal momentum fraction x, and 
integrated over all transverse momenta k\. 

As graphically illustrated in Fig. 1, the coupled evolution equations for the parton 
densities Qi of quarks, Qi of antiquarks, and G of gluons, respectively, can be summarized 
now as follows fill: 



| + |r) Qi(x,r,t) 



A + A'{G} + B'{Qi} 



Q, 



B + S qg {Q t } 



G 



i=i i=i 



(13) 



i + + B'{Qi} 



Qi + 



B + S.JQi} 



G 



+ E s g dQi} Qi + E 4(4} Qi 

3=1 3=1 

n f 

C + D + &{G} + J2{D'{QA + D'{Q 3 }} - S gg {G] 



(14) 
G 



3=1 



n f r 

+ E 

n / r 

+ E 

i=i 



E + E'{Q,} + S gq {G} 



E + E'{Q 3 } + S gq -{G} 



Qi 
Qi 



(15) 



Note that this set of equations describes the time evolution of the parton (longitudinal) 
momentum distributions Q = xq, Q = xq and G = xg, rather than of the parton num- 
ber densities q, q and g. The left hand side of these equations describes just the free 
streaming in the absence of interactions, whereas on the right hand side the various inte- 
gral operators acting on the parton momentum densities desribe the change of the densities 
due to branching, fusion, and scattering processes. In contrast to the branching operators 
(A, B, . . . , E), the operators that describe fusions (A'{Q}, B'{Q}, . . . , E'{G}) and scattering 
processes (S qq {Q}, . . . , S gg {G}) are functional of the nuclear parton distributions (labeled 



S 



with a ' ') of quarks Q (Q) and gluons G. 



In order to represent the integral operators in the evolution equations (|i3|) -([i5|) in a 
compact form, let me introduce effective coupling functions defined in terms of the running 
QCD coupling strength in one-loop order 

127T 

" s(g2) = (33-2n / )ln(Q 2 /A 2 ) ' (16) 

where A is the QCD renormalization scale and nj is the number of quark flavors that can 
be probed at scale Q 2 . The following definitions will be employed: 



OL* IT = 



Jo 1 dx (g{x, t, t) + J2j [Qj (x, t, t) + qj (x, r, t)]) 



i{r) = ^ , ((x,r) = 



(17) 

M 2 = \x\P/t 



2ttt ' ^ v ' ' M 2 

ivi~ — \j;\jr 1 7 

Here a s (r) describes the x-averaged QCD coupling at scale Q 2 = \x\P/r, and the age 
r represents the mean life-time of a parton with longitudinal momentum fraction x and 



virtuality Q . Finally I denote by "y a ^bc(z) the usual branching functions [|13|,|14 

~2\ 



l + z 2 



lg ^ gg {z) = 2 C G (,(1 - Z ) + Z J- + — £) 

l 9 ^{z) = \ (f + (1 - *f) , (18) 

where Cq = = 4/3, Cq = n c = 3, and z = Xb/x a is the fraction of x-values of daughter 
to mother partons in the branching a — > 6c. 

As is evident from the evolutions equations (|l3l) -(|T5|), the rate of time change of the num- 



ber of partons receives various contributions from the different interactions of the cascading 
partons, represented by their momentum densities Q, Q and G, from 1 — > 2 branchings, 
2 — >• 1 fusions and 2^2 scatterings. In addition to the well known parton evolution in 
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vacuum (e.g. jet development in e + e _ -annihilation), in which only spontanous branching 
processes contribute, in the present scenario of parton evolution in nuclear matter, the in- 
teractions with the medium, i.e. with the nuclear partons, give rise to fusion and scattering 
processes, and indirectly also to stimulated branchings. The physical picture underlying this 
probabilistic approach to generate such parton cascades from sequential elementary 1 — > 2, 
2 — ► 1 and 2^2 interactions is a direct generalization of the standard perturbative QCD 



parton shower picture |15| , |23| . For more details I refer to Ref. [[y]. Here I will only list the 
corresponding expressions of the integral operators in eqs. (|T^)-([T5|) corresponding to the 
diagrams in Fig. 1. 

2.2 Branching processes 

For cascading quarks the 1 — > 2 branching processes are balanced by the loss (— ) and gain 
(+) rates 



AQ l 



dz 



Qi(x,T,t) - Qi[-,T,t 



X 



X 



+ BG = + J o dzG{-,T,t) £(r) lg ^ qq {z) . 



(19) 



The corresponding branching rates for antiquarks are obtained upon substitution Qi <-> Qi. 
The loss and gain branching rates of gluons are given by: 



dz 



2 G(x,r,t) -G[-,r,t 



n f G(x, r, t) / dz f (r) J g ^ qq -(z) 



- CG = 

- DG = 

+ EQ j = + J dz Qj [^,r,t\ £(t) j q ^ gg {z) 

+ EQj = + J dz Qj (^> r >*) £( r ) 1q^qg( Z ) 



£( r ) lg^gg( z ) 



(20) 



2.3 Fusion processes 

The 2 —>■ 1 fusion rates that describe the effective loss of quarks due to fusions are: 



A'{G}Qi 



Pn c, 



gg~*g 



dz 



Qi{x,r,t)G (^^J C(- Z ,r 
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Qi (xz,r,t) G(x(l - z)) 



- B'iQ,} Qi = - PN c qq _ g £ dz Qi(x, r, t)Qi ^ - ^ ig-^z) C rj . (21) 

Here p^ is the nucleoli density of the nuclear medium (to be specified below), and the 
constants c qg ^ q = c qg ^ q = 1/8, c qq ^ g = 8/9 arise from the difference of flux factors (color 
and spin) for fusions compared to branchings [p5|| . Similarly the fusion rates of antiquarks 
are given by replacing Qi <-> Qi. The corresponding loss and gain fusion rates of gluons are: 

- C'{G}G = -p N c„ jf 1 dz[G(x,r,t)G {^^j C(f,r) 

- G (XZ, T, t) G (x(l - Z)) ( (X, T) lg->gg(z) 

* ~ ( x{\ — z)\ fx 

- D '{Qj} G = - p N c gq ^ q G(x, r, t) dz Qj ( J j q -+ gq (z) ( f -, r 

- D'-fQ,-} G = - p N c gq ^ q G(x, r, t) y ( J 7,_> OT (z) C ( -, t 

- 1 r 1 

+ ^'{Qj} <2j = + 2 ^ Cqq^g J Q dz Qj (xz) Q. (x(l - z)) l g -> qq (z) C (x, r) 

. ~ _ 1 A 

+ E'iQjjQj = + -p N c qq ^ g j dz Qj (xz) Qj (x(l - z)) -f g ^ qq (z) ( (x,r) . (22) 



The flux factors are here c gg ^ g = 1/8, c gq ^ q = c gq ^ q = 1/8, c qq ^ g = 8/9. 

In (pi ) and ( 22) the argument x of the nuclear parton distributions Q, Q and G is 



understood as \x\, since the nuclear partons move in — P direction [c.f. eq. ([!])], but the 
measured nucleon structure functions are usually defined only for positive calues of x. Also 
note that in fl2"2|) the second term on the right in the expression for C'{G} G would be 
multiplied by a factor 1/2 if G and G would be treated on equal footing, because nuclear 
and shower gluons would then be indistinguishable. 

2.4 Scattering processes 
The 2 — > 2 collision rates for elastic scatterings of the cascading partons carrying momentum 
fractions X\ and x^ with the partons in the nuclear background medium of nucleon density 
Pn carrying momentum fractions x[ and x' 2 are represented by the integral operator 

f° dx 2 t, N f , 2 da ab ^ ab B(x,T,t) 



S ab {A} B = - f ^ i(x 2 ) / dpi Pjv 
j-i xo j dv\ 



'-1 X2 J dp ± X 
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\x\P f° dx 2 r, x da ab ^ ab 
\ 5- / A(x 2 



t 2 J -i X2 dp\ 

B(x 1 ,n,t) 



2 \x\P 

P1 = U F- 



x p N / t/n 



1 



T Ixl 



Xi 

/° rfX2 At \ f XlP/T A 2 d °ab-,ab B(x 1} T,t) , 



where A,B = G, Qj, Qj and the labels as before the nuclear parton distributions, whereas 
the distributions without ' ~ ' refer to the cascading partons. As in the case of fusions, eqs. 
(f2l|) and (l22|) ) the value of x 2 in the argument of the nuclear parton distribution A is meant 



as I £2 1 when employing standard parametrizations of nucleon structure functions. 

The momentum fractions of two partons before (£1,2:2) and after (x' l ,x' 2 ) scattering 
involving a relative transverse momentum exchange p\ are related by 

x i = x 'i + 1 , P± ^ P 2 > x 2 = x' 2 + - - P± , (24) 

[x'i — x 2 )P l {x 2 — x x )P 



from which follows the value of X\ at which the function B(xi,r,t) in eq. ( 23 ) is to be 
evaluated, 

2 

X l = X + ( P \p 2 • ( 25 ) 

[x — X 2 )P Z 

For the quarks (and similar antiquarks) one has the processes q^g — > qig, q^qj — »■ qiqj 
and qi(jj — > q^j, whereas for gluons the contributing processes are gg — > gg, gqj — > gqj and 
gq~j — > gqj. The parton-parton cross-sections da a b->cd/ dp± in the scattering function (|23j) are 
related to the squared scattering amplitudes \Ai a b^cd\ 2 , averaged over initial spin and color 
states and summed over the final states, by (neglecting the quark masses), 
da ab ^ cd (s,pl) _ ™ 2 s (p1) |X t ,2 

^2 ~~ ^ab^cd p \Mab-*cd\ , l^Oj 

where the degeneracy factors V ab = (l + 5 a &) _1 accounts for the identical particle effect in the 
initial state if a and b are truly indistinguishable, and correspondingly V c d is the statistical 
factor for the final state. However, since cascading particles and nuclear partons are assumed 
to be distinguishable, one has V ab = 1 always (but not so for T> c d). Since needed later, let 



me list for completeness the relevant squared matrix elements for the various processes |26 
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The variables s, t, u are the kinematic invariants of the parton-parton scattering with 
s + t + u = 0, and p\ = tu/s for massless particles. For massive quarks the corresponding 
scattering matrix-elements can be found in Ref. . 



3. ANALYTICAL SOLUTIONS FOR THE PARTON SHOWER FUNCTIONS 



In the following I will present a method for solving the evolution equations (|l3|)-(|l"5|) 
analytically and will investigate some characteristic features of the solutions. Recall the 
scenario from Sec. 2 with a highly energetic incident quark or gluon, produced at some point 
of time to inside the target nucleus, with longitudinal momentum fraction xq ~ 1, transverse 
momentum k± = 0, and virtuality Ql ^> A, that initiates a cascade (or shower) of secondary 
quarks, antiquarks, and gluons. By solving the evolution equations, one can follow the time 
evolution of this parton shower in the nuclear background medium consisting of the partons 
of the heavy nucleus. In the present section I will solve the evolution equations (]T3|)-(|T5|) 
for the parton momentum densities Q (Q) and G, or equivalently the number densities q (q) 
and g. The synonym parton shower functions as generic term for these distributions will be 
frequently used. As I will show, the evolution equations can be solved in closed form, if one 
restricts to a somewhat idealized scenario that allows for a transparent analysis (without 
however giving up the physical relevance of application to realistic situations). 

3.1. Assumptions and approximations 

The assumptions and approximations that I will have to impose in order to achieve 
tractable solutions are the following: 
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(i) As before I will employ a rigorous distinction between the cascading partons produced by 

the evolving parton shower and the partons of the nucleus which represent the medium. 
Only the shower partons are followed dynamically. The nuclear parton distributions 
are approximated by a constant behaviour, motivated by a rather weak x-dependence 
of the nucleon structure functions at small values of x [ ffijfl . I will restrict the shower 
development to the small x-region where the probability for an interaction of a parton 
is the largest, because of the logarithmically large integration over the momenta of the 
involved partons and the large available phase-space. 

(ii) The nuclear background medium is assumed to be of very large (ideally of infinite) size 

with a uniform spatial distribution of nuclear partons, so that one can ignore boundary 
effects and average over the spatial coordinates. Also, I make the so-called "sudden 
approximation", i.e. the nuclear parton distributions remain essentially unaltered by 
the interactions with the shower partons, at least on time scales characteristic for the 
" intra-cascade activity" . However, once a nuclear parton has interacted with a shower 
parton, it is considered to be materialized from the coherent medium and becomes part 
of the evolving cascade. Thus the medium acts as a particle and energy-momentum 
reservoir. 

(iii) For the nuclear partons of the medium, I will ignore effects as nuclear shadowing, 
color coherence, long range correlations, etc., that affect the momentum and number 
distributions of those quanta. I approximate the nuclear distributions as simple su- 
perposition of the experimentally measured nucleon structure functions and consider 
the nuclear parton density as a function of the number of nucleons A per unit area. 
The contributions of valence quarks are neglected, that is, the matter is taken to be 
effectively baryon free, an assumption that is reasonable for the valence quark depleted 
small x-region. 

Accordingly, the attention is restricted here to the diffusion of cascading partons in 
longitudinal direction and the dissipation of longitudinal momentum (or energy) of the 
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shower particles. That is, I will solve for the time evolution of the parton shower functions 
as a function of longitudinal momentum fraction x and age r. A more realistic evolution 
including also the diffusion in transverse momentum k±, the lateral spread in fj_, as well as 
equal treatment of shower and nuclear partons, requires a numerical analysis, which will be 
addressed in a separate work. 

3.2. Simplified form of the evolution equations 

Due to the dependence of the integral operators for the fusion processes and the scatter- 
ing processes on the parton distributions Q, Q, and G of the nuclear background medium, 
the right hand side of the evolution equations involves quadratic forms of the parton dis- 
tributions, because of the coupling of the shower partons to the nuclear partons. A general 
analytical solution of this non-linear problem is therefore not possible. To overcome this 
obstacle, I explicitely distinguish these two sources [c.f. item (i) above]. I approximate the 
nuclear parton densities as scale {Q 2 ) independent and neglect their time variation, as well 
as nuclear shadowing effects, etc., [c.f. item (ii)] and simply represent them by the nucleon 
structure functions multiplied with a nuclear factor, 

Qi(x) = Pn xqi(\x\,Ql) , G{x) = p N xg(\x\,Ql) , (28) 

where g denote the usual (time-independent) measured nucleon structure functions, in 
principle to be evaluated at the initial scale Qq at which the primary parton is produced. 
It is important to realize that in ( p8[) the values of x = p z /P < in accordance with eq. 
(H), i.e. the fact that the nuclear partons move in — z-direction. Therefore, although the 
number densities cji, g are (as probability densities) necessarily always positive definite, the 
momentum densities are here Qi,G < 0. The nuclear factor depends on the geometry 
and the nucleon density p5| , p9| , 



3 

p N = -R A n N , (29) 

where Ra = r A 1//3 , r = 1.2 fm, and = A/ (4tiR\/3) are the nuclear radius and 
the nucleon density of the traversed nucleus, respectively. Furthermore, motivated by the 
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aforementioned weak x-dependence of the nucleon structure functions at small values of 
x<10~ 2 (see however pBfl), one may approximate the nuclear parton distributions (^) by 
a constant Q°, respectively G°, times the nuclear factor (|29|): 



q a 1 / 3 ~ q a 1 / 3 ~ 

Qi(x) = Q^x) = - — — Vl Q° , G(x) = - — — G° . (30) 



57T Tq 57T Tq 

Here the overall minus sign ensures that the nuclear momentum densities are negative def- 



inite, as explained after eq. (p28|). These x- and t-independent parametrizations mimic in a 
crude approximation the effective number of nuclear partons per unit area, weighted with 
their longitudinal momentum, as 'seen' by the evolving shower partons. Since at small x the 
contribution of the valence quarks is negligable, the distributions Q (Q) are understood to 
represent solely the sea quarks (sea antiquarks) for which I take rif = 3, that is, I account 
for the flavors i = 1, 2, 3 = u, d, s with the relative proportions 



Vi 



1 for i — u (u) 

1 for i = d{d) ■ (31) 
0.43 for i = s(s) 



As reasonable values for Q° and G° one may take: 

Q° w 0.2 , G° w 3 . (32) 

As long one is concerned with color and flavor insensitive kinetics of the parton evolution, it 
is sufficient and convenient to consider n/ = 3 identical types of quarks and antiquarks and 
the only flavor dependence comes from the relative quark admixtures ([H|) in the nuclear 
medium. Due to the absence of valence quarks the baryon number is zero and particle- 
antiparticle symmetry implies Qi = Qi. Therefore it is convenient to define the flavor 
singlet combination by the sum over all quark flavors, 

Q := E (Qi + Qi) ( 33 ) 
and consider quarks and corresponding antiquarks as one particle species. 
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Collecting all these ingredients, one can now rewrite the evolution equations (|13l)-(|T5|) in 
a much simpler, linear form: 



|Q(*,r,t) 



A + A'{G } + \b'{kQ°} - n f S qq {KQ } + |- 



Q 



+ 



2B + S M {/<Q } 



G 



— -Fgg Q + -Foo G 



</.</ 



(34) 



3_ 

at 



G(x, r, t) 



G + D + G'{G } + n/L^Q } - S gg {G } + 



d_ 



G 



+ 



E + -^{«:g } + S gq {G } 



Q 



(35) 



Here rji is defined by (|31|) and the constant k accounts for the effective number of quark 
flavors in the nuclear medium that are probed by the cascading gluons, 



— V* 



0.81 



(36) 



In Q3"3| ) and ( P5] ) both the parton densities and the integral operators depend still on the 
three variables t, r, and x, corresponding to the time development of the vituality (age) 
distribution and energy (longitudinal momentum) distribution, respectively. For each age 
t = \x\P/Q 2 one must in general expect a different x-dependence (and in principle also flavor 
dependence), reflected by an "age distribution" A a (r, t;x) which represents the probability 
that a parton a = q,g with x at time t has an age r, and which, for quarks and gluons, 
respectively, is defined through M 



Q(x,r,t) = A q {r,t;x)Q(x,t) 



with normalization J °° drA a (T, t; x) = 1, and 



A a (r,t) = / dxA a (r,t;x) 



G(x,r,t) = A g (r,t; x) G(x,t) , 



(37) 



fit) = dr dxA a (r,t;x) f(x,r,t) . (38) 
Jo Jo 

The age distribution evolves parallel with laboratory time in between interactions (free 
streaming), but due to scatterings it receives modifications because the age of the scattered 
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parton is reset to a younger age corresponding to a larger virtuality caused by the momentum 
transfer of the interaction. In the absence of a medium A = S(t — r). Therefore the age 
distribution reflects the influence of the nuclear medium and determines the characteristic 
age r. The age r must hence be related to the mean free time in between scatterings and 
consequently introduces an external scale depending on the density of the medium. The 
separation of the medium influence in the parton densities by introducing A, suggests to 
describe the parton evolution by the time rate of change with respect to the variable x only 
and interpret the age distribution as externally induced feedback to the shower functions 
via ( |37D and (|34"D , ([35]). Therefore let me rewrite the latter evolution equations in terms of 
Q(x, t), G(x,t) and A(r,t), defined in (|38~D , by substituting the factorized forms (|37|) : 

+ Ms+IM G + ^ G (39) 

Here the terms involving A must in reality determined selfconsistently by the dynamic in- 
terplay between the shower particles and the nuclear environment. It contains implicitely 
information about the nuclear structure which at this point has to be provided as phe- 
nomeno logical input. However, the actual form of A is irrelevant if one is not interested in 
the age distribution itself, but rather is concerned with the time evolution of the partons 
with respect to x and average r. In this case one may multiply both sides of the equations 
by A, integrate over r subject to the conditions fl38|), so that the explicit dependence on A 
drops out. Then one is left with eqs. (|39|), (|iOD, now without the terms A" 1 (d/dt + d/dr) A 
(taking A q ~ A g ), but one has to replace everywhere in the operators F ab 

POO 

t — ► fit) = v T (t)t = / drrA(T,t) , (41) 

Jo 

where 

+® < 1 w 

is the "velocity" by which the typical f changes with laboratory time t. The evolution 
equations can now be written in a compact form [[HJ: If one represents the time evolving 



parton distributions in a symbolic matrix form as U(t,t = 0) = exp[— iHt] with H the 
Hamiltonian of the system, then the evolving parton state will form a vector (Q, G) which 
satisfies the following Schrodinger type equation 



d 








( 


-F 


F \ 

1 qg 


(e\ 






= -iH 








dt 


\ G ) 






\ 


p 

1 91 


~ F 99 ) 





(43) 



As I will now show, this matrix equation, or equivalently the simplified evolution equations 
([39]) and (f40|) averaged over r according to (fill) , can be solved in closed form by diagonalizing 
the Hamiltonian H, subject to some chosen initial condition at t = t . 

3.3. Solutions for the parton evolution 

The integral operators in the matrix equation fl43|) are now forms of the average age f(t) 
and the longitudinal momentum fractions x = p z /P of the shower particles only. As stated 
before, the time-dependent f(t) loosely speaking reflects the structure of the nuclear medium 
as probed by the shower partons, and it reduces in free space to the trivial dependence f = t, 
i.e. v T = 1 in (^2|). Thus, one qualitatively expects a delayed parton evolution in medium 
as compared to vacuum, because repeated rejuvenations of the shower partons must lead to 
f < t, or v T < 1. If v T (t) has only a weak time-dependence, v T ~ const., then an appropriate 
solution ansatz for the parton shower functions is the following factorized form with a power 
dependence in 1/x and an exponential decay with some time-dependent function x(t) : 

Q(x,t) = J™ drA q {r,t) Q{x,r,t) = J\f q ( 6(x) + \x\ s 0(-x)J expf-^xW] 
G(x,t) = J™ drA g {r,t) G(x,r,t) = M g ( \x\~ s 6(x) + \x\ s 9{-x)j exp[-px(t)] ■ ( 44 ) 

Let me explain the physical motivation of this ansatz: 

(i) The introduced power s is understood to be positive definite, s > 0, and the constants 
M q and M g are normalization factors. As will become clear, the 'parameter' s plays the 
role of a conjugate variable to the variable x, controlling the longitudinal momentum 
distribution in x. As the shower evolves, s will change with time, thereby correlating 
the x and the t dependence of the shower functions. 
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(ii) The form of the ^-dependence in brackets is motivated by the particular dynamics 

viewed in the CM^n frame (Sec. 2.1) in which one naturally expects some "drag 
effect" of the (in — z direction moving) nuclear partons on the evolving parton shower. 
Although the shower partons initially propagate in +z direction, one must account at 
some point for particles moving in both +z and —z direction. Shower partons moving 
in +z direction will progressively degrade their longitudinal momentum by branchings 
and scatterings, and will eventually reverse their longitudinal direction, dragged along 
and accelerated in —z direction by absorption of (fusions with) comoving nuclear 
partons. 

(iii) The exponential factor accounts for the fact that the running coupling strength a s in 
the operators F a b of fl43|) depends on f(t) and therefore implicitely on t. By introducing 
the function 

x(t ) S /V^ = , (45) 

J t 27rr(r) 2nbv T \a s {v T t ) J 

where b = (lln c — 2n/)/127r = 0.72 and f = v T t with v T assumed to be only slowly 
varying with t, it is possible to absorb the medium dependent factor a s (f) / (2nf) in the 
exponential time dependence of the distributions (f|4|). Note that for fixed coupling and 
v T = const., one has x °^ hi^Ao)- Even in the general case x represents approximately 
the logarithm of the laboratory time t. The real time dependence in the laboratory 
frame is regained by the inverse relation 

xP 



cxp 



In (ff) exp(-2irbv T x) 



,v T A 2 , 

where Ql is the virtuality of the shower initiating parton at t . 



(46) 



Substituting the ansatz ([0]) into the time-evolution equation (^3[), one arrives at the follow- 
ing eigenvalue equations with respect to the new time evolution variable x{t)i corresponding 
to the diagonalized Hamiltonian H: 
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d 








dx 


\ G j 




\ G ) 



(47) 



or equivalently 



li = F qq (s 



- a V = - a Fgg(s) ~ F gq ( S ) 



(48) 



Here the scalar functions F a b(s) correspond to the operators in the Hamiltonian, F a b, defined 
by © and 



1 



F K (s) = A(s) + A' m (s) + ^B' {kQ0} (s) - n f S\^ 



2 {«QT 



F gg (s) = 2B(s) + S[l% } (s) 

F gg (s) = n,E + n f E' {KQ0} (s) + n f S$ } (s) 

F gg (s) = C(s) + C' m (s) + D(s) + n f D' {K Q 0} (s) - S<gJ } ( a ) . (49) 



arc 



The branching kernels in 

A(s) = £dz (l-^) 7«-h»(*) = C Q ^2ij(s + 2) + 2 lE -j-^ { 
B(s) 



l)(s + 2) 



3 



1 , if s 2 + 3s + 4 

o l9 ^ [Z) " 2 \( a + l)(a + 2)(s + 3) 



C(a) 



= C G {2^(s + 2) + 2 7£ - 4 



s 2 + 3s + 3 



E{s) 



Cr 



s(s + l)(s + 2)(s + 3) 6 
s 2 + 3s + 4 



11 + n f 



3C, 



G 



(50) 



io - s(s + l)(s + 2) ' 

where the branching functions 7 (W & C (< 2 ) are given by flTB]), and as before rij = 3, Cq — ,, . 



ni-l 



Cg = n c with n c = 3. The Psi (Digamma) function is defined as ip(z) = d[lnT(z)]/dz and 
je — 0.5772 is the Euler constant. Note that these branching kernels are finite, i.e. free of 
infrared singularities, because the infrared divergent contributions cancel explicitely due to 
the balance of gain and loss terms. 

Similarly the fusion kernels are evaluated as: 
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A \g»}( s ) = - G° h A (v T )c 



dz (z- z s ^j %-> qg (z) 
G° h A (v T ) c„ C Q I 2 if,(a + 2) + 2 lE - (g - ^ - 2) - ^ 



S {reQ°}( S 2 



- KQ h (V T ) Cgg^g 

CU Js) + n f D' d0 Js) = -& h A (v T ) c 



dzZ ^g^ q q(z) = - kQ° h, A (v T ) C 



qq^g 



■^99^9 



dz (z-Z^j lg^ aa {z) 



- kQ° h A {y T ) Cgq—,q n f f dzz^ q ^ gq {z) 

J 

& h*M o^ 3 C G { 2 + 2) + 2 7E - 4 , (<+1)( , + 2)(< + 3) 



s 2 + 3s + 3 



11 /^c^4C7q^ ■ 
6 G° c gg _ 5 3C G n/ 

K<5° h A {v T ) Cgg^g 



- kQ° h A (v T ) C 



qq^g 



dz Z S ^g^qq(z) 

1 s 2 + 3s + 4 

2 (s + l)(s + 2)(s + 3) 



(51) 



The constants c a b_> c in ([51]) are given after eqs. ([21] ) and ([22]). The dimensionless function 
/i' 4 contains the nuclear dependence, eqs. fl29|) and (|30|), 



|x|P 



r§A 2 



exp 



- In ( J exp(-2vr6t; r x) 



(52) 



Finally the scattering kernels are obtained as: 



7(99) 



1 



- « s (V T t) KQ h A (v T ) T.qq^qq 



= - a s { Vr t) G° h A { Vr ) ^gq 



qg^qg 



S ( g (s) = - a s (v T t) G° h A (v T ) S„ . (53) 



{G }^ 4 

These scattering kernels result from the operators (^) in the limit of the "small-angle 
approximation" or "Landau approximation" [ 30| , in which one restricts to scatterings with 
small transverse momentum exchange pi <C s < P 2 . In this case the scattering operator 
( p3|) simplifies considerably to 



S ab {A}B 



00 r 1 dxi r° dx 2 f 00 , o n/ s 7/ \ 
«Ti / / / dp ± B(x 1 ,T,t) A(x 2 ) 



x 



Xi J -I x 2 Jo 

d&ab^ab 



„2 



5 (a/ 2 — x) S 



\x\P s 

pi , 



(54) 
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where scattering at small angles implies [c.f. eq. (^4 



x 1 = x [ + — 



Pi 



X 



x 2 = x' 2 + — 



A 



X. 



2 5 



(55) 



{x' l -x 2 )P 2 1 ' 1 (x' 2 - Xl )P 2 

and the parton-parton cross-sections da a b^ab/ dp\ given by and fl2"T|) can be replaced by 
their small-angle contributions 



da, 



ab—*ab 



dp] 



no; 
P 4 ± 



-'ab—tab 



(56) 



where the constants S a fc^afe follow from the squared matrix elements \A4 a b-, a b\ 2 , eq. 
by observing that at small angles p\ = — i <C s and keeping only the leading terms in p\/s: 



y = y 

^qq^qq — ••-./../, •</.</; 



y - y~ - - - 

^qiQi^qiqi ~ ^qiq^qiqi 



J qiqi^qiqi 



8 
9 



y = y - y~ - -9 

^qg^qg — ^qig^q t g ~ ^qig^qtg ~ * 



y = y - y~ - 9 

^gq-^gq — ^gqi^gqi ~ ^gq^gn ~ * 



J gg->gg 



(57) 



Returning to the eigenvalue equations (58) and combining them, one obtains 



H + A+(s) 



H + A_(s) 



which implies 

A + (s) = -i \f w (s) + F gg (s) 
A-(-s) = — F w (s) + F 99 (s) 



[i^W - Fgg(s)Y + 4F qg (s)F gq ( S ) 



and, corresponding to these functions, the constraints 



M q F gg (s) + A H 



for 



or 



Fgg(s) 



M q F gg (s) + A_(s) 
Thus the solutions are given by 



for fi 



-A_ . 



(58) 



(59) 



(60) 



(61) 
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Q(x,t)=Af q (£) ex V [X + (s) X (t)} 
M q F gq (s) fl- s 



« = W^1;J exp[A +(S )x(t)] , (62) 



or 



Q(x,t) = A/" 9 Q exp[A_(s)*(t)] 

C ^"- ^TTM (§)' "PlM')*ffll ■ (63, 

or by a linear combination of these two solutions. In general one can write the solutions for 
the parton number distributions p = Qj 'x = J2i (<?i + Qi) and g = G/x as: 

dx 



p(x, t)dx 



x s+l 



, \ j doc 
g{x, t)dx 



x s+l 



a e A+(»)x(t) + a _ e M«)x(*) 



Q+^ot e A+(») x (t) + Q - F gq e A_(s)xft) 



(64) 



Fgg + A+ Fgg + A. 

Fig. 2 shows the behaviour of the functions F qq , F qg , F gq , F gg and Fig. 3 of A+, A_, and 
the derivatives A^, A" with respect to the variable s. The three CcLSCS cl)-C ) correspond to 
different nuclear density of the medium as controlled by the nuclear mass number entering 
the function h A , eq. fl52|): a) A = (free space), b) A = 50 (moderately dense medium), 
and c) A = 200 (dense medium). Evidently A + vanishes at s = Sq ~ 1.0; 1.1; 1.2 for the 
cases a); b); c), respectively, and is positive for s < so, and negative for s > sq. From ( |64"D 
it is obvious that the number distributions of shower partons behave as dx/x s+1 , so that 
one can interpret the region with s < sq as the developing stage of the shower, with s = sq 
representing the situation at the shower maximum, and the region with s > Sq as the fading 
stage leading to the tail end of the shower. The variable s therefore indicates the shower age 
at a given value of x: in a 'young' shower (s < sq) the number of shower particles increases 
with age, while in an 'old' shower (s > sq), it decreases. 

Solutions corresponding to a certain initial condition are obtained by summing up these 
parton shower functions with respect to the variable s, weighting them with an appropriate 
amplitude. This procedure is known as the Mellin transformation method. It has been 
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applied extensively to similar problems of electron-photon showers in the early days of 
cosmic ray theory |31|j32]j • The reason for introducing Mellin transforms of the parton shower 
functions is that they are more easily determined than the original functions. Let me denote 
by M q (M g ) the Mellin transform of the function p (g), 

MJs,t) = / dxx s p{x 1 t) , MJs,t) = / dxx s g(x,t) , (65) 
Jo Jo 

then the inverse transformation is 

1 rc+ico MJs.t) . , 1 /- c + i0 ° , M Q (s,t) 

where s is now a complex parameter, and the path of the integral is running paralell to the 



imaginary axis. By analytic continuation of the shower functions ( 64]) to complex values of 
s, one obtains the solutions 



rc+ioo | 

fpfat) = t^—. ds- 



2iri Jc-ioo x s+1 

I rc+ioo | 

2m Jc-ioo x s+L 



a e A+(*)x(t) + a _ e*-Wx(t) 



(67) 



where a + and a_ are functions of the complex parameter s, and are determined by the initial 
conditions imposed on p and g. 

I will now discuss the two cases of 1) a parton shower initiated by a quark, and 2) a 
parton shower triggered by a primary gluon. In accord with the physical picture that was 
outlined in Sec. 2, the shower axis z is defined along the momentum of the primary quark 
or gluon, so that the incident parton has zero transverse momentum. Since I consider here 
the longitudinal evolution only, the time development of the parton cascade proceeds by 
increasingly populating the region with small values of x as the shower develops along the 
shower axis which corresponds to growing values of s. According to eq. @, the initiating 
primary parton carries a longitudinal momentum p z0 = XqP and has a initial time-like 



virtuality Ql such that its energy is determined by E = J (xqP) 2 + Qq = XqP + QI/(2xqP), 
where Q <S P is assumed. 
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3.3.1 Parton shower initiated by an incident quark 

If one considers a parton shower initiated by a primary quark, then the initial conditions 
at time t — to are 



p(x,t ) = 5(x Q -x) 6 [E - -Py^o + pf 



Inserting this initial condition into fl65|), one has at to = 0, 

M q {s,t ) = x s , M„(s,t n ) = 

which means that 



sOMo) 



. 



(68) 



(69) 



] rc+ic 



x / x 



g(x,t = 0) = . 



(70) 



The parameters a+ and a_ are now completely determined by (|68|)-(|70|), together with 



a + + a_ 



x 



a+F ( 



o j 



+ 



Fgg + ^+ Fqo + A_ 



. 



(71) 
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If one introduces functions H + and H_ through 

a + = H + (s) Xq , a_ = H-(s) x\ 



) 



(72) 



or 



H 4 



F gg (s) + X + (s) 



F gg (s) + A_(s) 



(73) 



X+(s) - A_(s) ' v y A+(a) - A_(a) ' 

one can now express the solutions for the differential spectra of secondary quarks and gluons 
produced by a primary quark as 



p(x, t)dx 



1 


dx 


rc+ioo 


ds 


(Xq 


2iii 


X 


J c—ioo 


K x 


1 


dx 


rc+ioo 


ds 


(Xq 


2m 


X 


J c—ioo 


V x 



H+(s) e x+is)xit) + H(s)e x - is)x(t) 



e M<0x(t) _ e M*)xW 



(74) 



where 



L( S ) 



y/sF gq (s) 
X + (s) - A_( S ) 



(75) 
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Correspondingly, one can evaluate the integral spectra of secondary quarks and gluons, that 
is the integrated number densities of partons with longitudinal momentum fractions greater 
than a given value x, 

1 r c + io ° ds fx 



Ux,t) = / dx'q(x',pl,t) = — — ) H + (s)e x ^^ + HJs) e A -«x(t) 

Jx ~ 2m' Jc-ioc s \ x J L 



rxo | r 

I g (x,t) = / dx' g(x',p 2 ± ,t) = — 

Jx ~ 2m Jc 



1 fc+ioc ds fX \ 

— L ( S ) 



s 3 / 2 \ X 



e x+(s)x(t) _ e A-00xM 



■ (76) 



3.3.2 Parton shower initiated by an incident gluon 

If one considers instead of a primary quark an incident gluon, the initial conditions are 



p(x,t ) = , g(x,t ) = 5(x -x) 5 \E - pJxl + Qt 



(77) 



the procedure is completely analogous. In this case one has instead of (|74|) for the differential 
spectra: 

1 dx rc+i°o 



, , , I dx f c ^°° , fx y r— . 
p[x, t)dx = / ds — ysJ{ 

2m X Jc-ioo \ X J 



,M«)x(t) 



, . 1 dx f 

g{x,t)dx = — — / 
2m x Jc 



ds { — 

x 



jj_(s) e M«)x(*) + HJs) e x - {s)x{t) 



(78) 



where 



sXJs) - XJs) 



J(s) 



For the integral spectra one has instead of (ff6[) : 



(79) 



Jq(x,t) 

J a (ar,t) 



c + io ° ds (xq 



2m . 'c— zoo \/ S \ X 

1 r c + io ° /a; x " 

27U Jc-ioo S \ X 



J(s) 



3 A+( S )xW 



3 A-( S )xW 



#_( s ) e A+Wx(t) + # + ( s ) e A_(-)x(*) 



(80) 



3.3.3 Saddle point evaluation of the parton spectra 

The integrals F = p,g, I q , I g can be evaluated by numerical computation of the contour 
integral in the s plane, but alternatively one can also obtain analytical forms by employing 
the saddle point method which I will outline now. First note that Fig. 3 exhibits the fact 
that A + > A_ for all values of s, so that one can neglect the second term in these integrals 
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(oc exp[A_x]) compared to the first term (oc exp[A + x]). Therefore, dropping the second 
term the typical structure of the integrals is of the form 

F ( x ,t) = — ( C+lC ° dss- k (-Yn(s) e^+WxW 

] rc+ioo 

= — / dsTC{s) exp[X + (s)x + ys - kins] , (81) 

where the symbol TC under the integral stands for the functions H + , H_, L, J, the variable 
y = \n(xo/x), and the values of the power k are -1/2, 0, 1/2, 1, or 3/2, as can be read off 
the formulae of Sees. 3.3.1 and 3.3.2. Now, the functions TC are slowly varying with s, so 
that the integrand in (|8~ID is a product of a function TC with a weak s-dependence times an 
exponential of the form exp[0(s)], where 

<f)(s) = \+(s)x + ys - kins . (82) 

If d(f)/ds = and d 2 <j)/ds 2 > 0, at some point s — s, then the integrand in ( |8"ID may be 
approximated by a Gaussian function on the path parallel to the imaginary axis crossing the 
real axis at s = s, with the slowly varying functions TC(s) treated as constants and evaluated 
at s. Thus, by expanding 0(s) around s = s where d(j)/ds = 0, one can evaluate the integrals 
by the saddle point method (see e.g. 



X J 



exp[X + (s) X ] , (83) 



with the saddle point s defined by the equation 



^ = \' + (s)t + y-* = , (84) 
OS s 

which establishes a correlation between y = ln(xo/x) and x> i- e - between the momentum 
distribution and the time evolution variable \. Applying this formalism, one obtains the 
following results for the parton spectra F — p,g, I q , I g (dropping the bar on the s, i.e. s = s 
is the saddle point): 

1) For a primary quark of longitudinal momentum fraction xq, 
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2) For a primary gluon of longitudinal momentum fraction Xq, 
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Recall that the time evolution variable defined in (pE5[), relates these spectra to real 

time t via eq. (pESP, so that t can be resubstituted if desired, provided v T in (fj2"|) is known. I 
would like to remind the reader that one of the advantages for using x , ec i- (H§)> as evolution 
variable rather than the laboratory time i, was to absorb the nuclear medium dependence 
and thereby obtain a relatively model independent description. Thus, notwithstanding the 
uncertainty of medium effects, the properties of the parton shower evolution are essentially 
controlled by perturbative QCD, but the connection to real time resides in the relation fl46[). 
From this point of view, the parton spectra ( j35|) and (|86"D represent the elementary solutions 
to the evolution equation (|43|) under the approximations made in Sec. 3.1 and summarize 
the longitudinal parton cascade evolution in nuclear matter within the present framework. 
It is evident that the variation of the number of shower particles with time follows an expo- 
nential law with exponent X+(s)x{t) with the (implicitely medium dependent) 'absorption 
coefficient' (— A + ) describing the variation with time and the power (— s) controlling the 
momentum distribution in xq/x. 
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(86) 



It is worth noting that the distributions p and g reduce to the well known approximate 
solutions to the Q 2 -evolution according to the DGLAP equations |18[ , if one sets the nuclear 
matter density equal to zero. Then r = v T t = t or v T = 1, because in the absence of 
scattering and fusion processes the particles can only evolve by successive branchings with 
the age r being equal to laboratory time t. The connection between time- and Q 2 -evolution, 
as discussed in Sec. 2.1, is given by the average relation @, Q 2 = xP/t. Then, considering 
the evolution of the gluon distribution alone, at small x the behavior is controlled by the 
branching kernel F gg (s) = C(s) near s = [c.f. eq. (pOf)]. Hence, taking only the dominant 
term in F gg into account, one has 
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where Co = n c as before. Then the saddle point evaluation yields (setting x = 1): 
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where the time variable Xi when re-expressed in terms of Q 2 rather than t, is given by 
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(89) 



2tt6 \a s (Ql) J 2iTb \\n(Q 2 /A 2 ) J ' 

with Ql = P/t Q and Q 2 = xP/t (t > t ). Eq. (||) together with (|89| ) coincides with leading 
logarithmic behavior of the gluon distribution within the LLA ||18||. 



Let me return now to the parton spectra (|85D and (^6|), and examplarily exhibit some 
characteristic features of the longitudinal parton shower evolution in nuclear matter, as 
reflected by the Xi x an d nuclear (A) dependence of the spectra. In what follows I chose 
the virtuality of the cascade initiating parton Q = 30 GeV. An increase (decrease) of 
Qo generally results in a larger (smaller) overall multiplicity of secondaries and also shifts 
the evolution toward earlier (later) times. The initial longitudinal momentum fraction was 
set to xq = 1, since the shower development depends only on the ratio xq/x. Finally the 
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QCD scale A that enters the running coupling a s was taken equal to 230 MeV. In order to 
exhibit the nuclear medium influence on the shower development, as controlled by eq. (|52|), I 
considered as before the three ) A = (absence of medium), b) A = 50 (moderately 

dense medium), and, c) A = 200 (dense medium). 

In Figs. 4 and 5 the parton number densities p(x,t) and g(x,t) from fl35|) and (|S6l), 
respectively, are plotted for x = 10 _1 and x = 10 -3 versus % (it is helpful to keep in mind 
that for slowly varying a s , one has approximately x oc hit). The essential observations are 
basically the same for both quark and gluon initiated showers: 

(i) For a given value of the ratio x, the change with "time" \ of the spectra p and g 
reflects the aging of the shower. The multiplicity initially increases (young shower), goes 
through a peak (shower maximum), and then decreases again (old shower). The time it takes 
for a parton shower to evolve from Xo = 1 downwards to some value x naturally increases 
with decreasing x and at the same time the average multiplicity of produced particles grows 
rapidly. (The integral spectra I q and I g , which are not shown here, exhibit an analogous 
behaviour.) 

(ii) Comparing the ) to c), the influence of the nuclear medium manifests itself in a 
slowing down of the shower development, particularly at later times (% >, 2). In other words, 
the aging of the cascading partons is more delayed, the denser the medium is. Related 
to that is the clear increase of the parton multiplicity at a given x with larger nuclear 
density, because repeated energy-momentum transfers by scatterings and fusions rejuvenate 
the shower partons and amplify the bremsstrahlung processes. 

In order to exhibit the x-dependence of the parton spectra at fixed value of x let me define 
the so-called normal spectra |31]] at the shower maximum. The normal spectra correspond 
to an optimum value x = Xo, which coincides with the value of x that makes the function 
exp[0(s)] of eq. fl52"| ) a maximum, because the other terms change slowly with x- Thus, Xo 
is defined by the equation 

d 

— exp[0(s)] =0 at x = Xo (90) 
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or 

a+00 + (y + KrWx - 7) (|) = ■ (si) 

Recalling that eq. fl53|) is satisfied, one has 

A+(s) = at* = Xo • (92) 

Hence, at s = so, 



y — k 

Xo = 



h, ( ^ ) /.■ 



(93) 



a; ( So ) a; ( So ) 

The normal spectra of p, g, I q , I g and the corresponding optimum values xo are easily 
obtained from eqs. fl85|) and (p6|) by putting s = s , A + = 0, and x = Xo- The values of s , 



given by A + (s ) = 0, and of A' + (s ) can be read off Fig. 2. 

The top part of Fig. 6 displays the typical shape of the normal spectrum around shower 
maximum, corresponding to s = So and x = Xo- It is interesting that, when plotted as 
a normalized probability distribution p(xa), this normal spectrum is unique: first, it is 
the same for quarks and gluons, second, it is independent of the type of primary particle, 
and, third, it does not depend on the nuclear medium density. On the other hand, the 
corresponding optimum values Xo vary significantly for the different cases, as is evident from 
the lower part of Fig. 6, which shows xo versus x for the three aforementioned cases a) 
A = 0, b) A = 50 and c) A = 200. This behaviour is to be expected, since Xo marks the 
characteristic time that it takes for the parton shower to build up the normal spectrum. 
Looking at Fig. 6, one can conclude: 

(i) With decreasing x the values of Xo increase, roughly as ln(l/x). The slope is the same 
for both quarks and gluons, but Xo is always smaller than xl, reflecting the fact that the 
gluons evolve faster due to their larger interaction and radiation probability. 

(ii) The effect of the nuclear medium is again a substantial dilation of the shower devel- 
opment, and thus toward the establishment of the normal spectrum. This is clearly reflected 
by the steepening of the slope of the time scale Xo when proceeding from the free space case 
a) to the dense matter case c). 
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5. SUMMARY 

In this paper I presented an analytical method to solve the QCD evolution equations de- 
rived in Ref. [ffj which describe the cascading of quarks and gluons in a nuclear environment. 
These evolution equations extend the well known DGLAP equations for jet development in 
hard QCD processes, to the time evolution of parton cascades in a nuclear medium, includ- 
ing in addition to spontaneous branchings also fusion, scattering and stimulated branching 
processes. In order to obtain tractable closed solutions to the evolution equations, a number 
of idealizing approximations were necessary. Therefore, at this point, this approach does of 
course not provide quantitative predictions of parton transport phenomena in a dense mat- 
ter environment. The main uncertainty is here the lack of detailed knowledge of the form of 
the medium influence on the cascade evolution. Since frequency and type of interactions of 
the partons with the medium must depend crucially on the nuclear density and structure, 
one must supply some external input. In order to minimize the amount of phenomenological 
input, the medium influence was therefore hidden into a time evolution variable xW> the 
relation of which with the laboratory time t however is determined by the specific properties 
of the nuclear medium. As advantagous trade-off one gains that the calculated properties 
and spectra of the cascading partons are a rather firm consequence of perturbative QCD 
proceses in a nuclear medium being characterized by a specific time variable \. 

Of prime interest in this work was the multiplicative behaviour of cascading partons at 
small x, since such soft particles each take away only a negligable portion of the total energy, 
but form at the same time the bulk of multiplicity. The characteristic properties of parton 
shower evolution inside nuclear matter with respect to the variables time (x), longitudinal 
momentum fraction (x) and nuclear density (A) can be summarized as follows: 

(1) With respect to the rate of time change at fixed value of x, a typical parton shower 
can be divided roughly in three stages: (i) young shower, during the early stage the 
shower particles rapidly multiply due to branching processes (mainly gluon emission); 
(ii) shower maximum: for a short time the number of particles in the cascade remains 
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almost constant; (iii) old shower, the increasing dissipation of energy due to branchings 
as well as a growing number of fusion and scattering processes reduce the number of 
shower particles at given x during later times. 

(2) The mean of the ^-distribution of shower partons shifts with progressing time to smaller 
values of x, but this development is significantly delayed in a dense matter environment 
as compared to free space. The x-spectra exhibit a unique form around the shower 
maximum, independent of type of parton and of the nuclear density. This normal 
spectrum is established after a certain optimum time X = Xo-, which measures the 
characteristic time scale that depends on the density of the surrounding medium: the 
denser the matter, the longer it takes for the partons to build up to their normal 
spectrum. Also, gluons always reach this point earlier than quarks, because of their 
larger cross-section and radiation probability. 

As said before, the approach presented here has no truly predictive relevance for experi- 
mental observables yet. Rather than that it may be a first step towards a more fundamental 
understanding of " QCD in medium" , using well developed methods of renormalization group 
improved perturbation theory for " QCD in vacuum" . I believe that the most important is- 
sues to be addressed in the near future within such a program are the following: 

(1) The rigorous distinction between the shower particles and the nuclear partons needs to 
be dropped. Instead both should be treated on the same footing in order to describe a 
self-consistent dynamical multi-parton system. This would then incorporate also time- 
and x-dependent nuclear parton distributions that receive modifications due to the 
response of the shower evolution. An interesting aspect of such a dynamical interplay 
between scatterings, multiplication and recombination of partons is the question of a 
possible saturation |33j of the parton number densities that may render a finite parton 



density as x — > 0. 

(2) In the context of nuclear structure, the relevance and the form of the age distribution A 
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of vitual shower partons needs to be studied, because it provides an interface between 
the nuclear medium and the response of the cascading partons to this environment. 
Therefore A may be experimentally accessible by comparing inclusive quantities asso- 
ciated with parton evolution in e.g. deep inelastic scattering experiments on hadrons 
versus heavy nuclei. This issue is of great interest from both theoretical as well as 
experimental side. 

(3) Interference effects which are important particularly in the small x-domain need to be 

taken into account systematically. These include two well known phenomena, namely, 
a) the coherence of successive soft gluon emissions in the parton cascade leading to 
destructive interference and a suppression of gluon radiation, and, b) the Landau- 
Pomerantchuk-Migdal effect, i.e. the interference of induced gluon emission amplitudes 
in multiple parton collisions, which also causes a depletion in the gluon radiation 
spectrum. 

(4) A possible extension to include the time development of the transverse momentum dis- 

tribution of partons as well as the lateral spread of shower particles would be desirable. 
This would then provide a solution to the cascade development in full four dimensional 
momentum space, in virtuality Q 2 , longitudinal and transverse momentum components 
x and kj_, respectively. 

(5) Parallel to the previous items a comparison with a full six-dimensional phase-space 

analysis, as e.g. possible with Monte Carlo simulations [Q], is of great importance 
to test the accuracy of the analytical solutions. Such a twofold investigation can be 
explorative for both analytical and numerical solutions of the evolution equations, 
in the sense that mutual feedback may lead to an improved understanding of the 
microscopic parton dynamics in nuclear matter. 
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FIGURE CAPTIONS 



Figure 1: Diagrammatic representation of the coupled evolution equations ([i~3|)-(|l5f) for a) 
quarks (and similarly antiquarks) and b) gluons. In contrast to the branching processes, 
the fusion and scattering processes involve interactions of a cascade parton with a nuclear 
parton (marked by thick lines). In a) the various diagrams correspond to the integral 
operators —AQi, BG, —A'Qi, —B'Qi, J2j[S qq Qj + S qq Qj], S qg G, and in b) to —CG, —DG, 
Ej[EQj + EQj], -C'G, -D'G, EA E 'Q 3 + E'Qj], E,[4A + 4<AL S gg G. 

Figure 2: Form of the functions F qq , F qg , F gq , F gg , eq. ([4~9|), with respect to the variable 
s. The three plots a)-c) correspond to different nuclear density of the medium measured in 
terms of nuclear mass number: A = (free space), b) A = 50 (moderately dense medium), 
and c) A = 200 (dense medium). 

Figure 3: Functional behavior of A + , A_, eq. (|5^|), and of the derivatives A' + , A" versus the 
variable s. The cases a)-c) correspond to Fig. 2. 

Figure 4: Primary quark. Time-evolution of the differential spectra for the parton number 



densities p(x,t) and g(x,t) with respect to the variable x{p)i according to (|85|), at fixed 
values x = 10 _1 and x = 10~ 3 . The virtuality of the cascade initiating quark was chosen 
Q = 30 GeV. Cases a)-c) as in Fig. 2. 

Figure 5: Primary gluon: Development with time x(t) of the parton number densities 
p(x,t) and g(x,t) for the case of a shower initiating gluon, eq. (jS6l), with initial virtuality 
Qo = 30 GeV and at fixed values x = 10 _1 and x = 10~ 3 . Cases a)-c) as in Fig. 2. 

Figure 6: Top: Unique shape of the normal spectrum around shower maximum in terms 
of the normalized probability distribution p(xo) — a ( x ,Xo)/ 1 dxa(x,Xo) where a = p,g. 
Middle and bottom: The optimum values Xo for gluons and quarks versus x, corresponding 
to the point of time x{t) when the normal spectrum at a certain value of x is established. 
The three sets of curves a) and c) refer to different nuclear density, as in Fig. 2. The 
virtuality of the primary quark, respectively gluon, was chosen Q = 30 GeV. 
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